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Abstract. This work presents a preliminary study of 147 relativistic mean-field (RMF) hadronic 
models used in the literature, regarding their behavior in the nuclear matter regime. We analyze 
here different kinds of such models, namely: (i) linear models, (ii) nonlinear a 3 + a 4 models, 
(iii) (7 3 + a 4 + (O 4 models, (iv) models containing mixing terms in the fields a and (0, (v) density 
dependent models, and (vi) point-coupling ones. In the finite range models, the attractive (repulsive) 
interaction is described in the Lagrangian density by the a (to) field. The isospin dependence 
of the interaction is modeled by the p meson field. We submit these sets of RMF models to 
eleven macroscopic (experimental and empirical) constraints, used in a recent study in which 240 
Skyrme parametrizations were analyzed. Such constraints cover a wide range of properties related 
to symmetric nuclear matter (SNM), pure neutron matter (PNM), and both SNM and PNM. 
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Quantum Hadrodynamics (QHD) is an approach based on quantum field theory much 
used in the description of hadronic environments, such as nuclear and neutron matter. 
It is based on local Lagrangian densities whose free parameters are adjusted in order to 
reproduce basic nuclear matter bulk properties at zero temperature. In general, nuclear 
matter is well described by different versions of the relativistic mean-field (RMF) models 
constructed via the QHD approach. 

In order to select a set of RMF models which better describe nuclear matter properties, 
as well as the physics of pure neutron matter, we submit 147 relativistic parameteriza- 
tions to 1 1 constraints, out of which 4 are related to symmetric nuclear matter, 2 to pure 
neutron matter, and 5 related to the symmetry energy that involve both symmetric and 
pure neutron matter. All these constraints were taken at the saturation point (po), except 
the constraints related with the band region. For more details about the constraints and 
the criteria used for approval see Ref. [1]. 

The relativistic hadronic models [2, 3, 4] used here are described by the following 
Lagrangian densities: 



(a) Nonlinear finite range (containing 123 models): 



where 
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(b) Density dependent models (containing 6 models): 
£dd = W(ifd^~M)Y + r (J (p)o\j/Y-T (D (p)\j7y^co^Y- 1J ^¥fPti-W 
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where r,-(p) = n(p )/;(x), with /}(*) = g f ' ' for / = (7, co, and 

1 + C/(JC + a/) 

r p (p) = r p (p ) e - fl ( x - 1 ) with^ = p/p . 

In both Lagrangian densities, M, m ; , with i = a, ft) and p are the nucleon and the 
mesons masses, respectively. The tensor fields are given by F^ v = dyCfy — d^COv and 

B^iv = dyp^ — dfxPv. gi (i = a, ft), and p), A, and B are the coupling constants, 
(c) Point-coupling models (containing 18 models): 

^pc = Wr%-M)^-^(w) 2 -y(W) 3 -f (w) 4 

- f (W» 2 - f (W» 4 - ^(Wfw) 2 - ^fiwfwf 
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where as, As, 75, ay, yy, ajy, Ytv, (Xts, ill, f]2, f]3 are the coupling constants. The 
subscripts identify the coupling: S stands for scalar, V for vector, and T for isovector. 

The equations of state for these groups, e.g., the energy density and pressure, are 
calculated from the energy-momentum tensor (in the mean-field approximation): 
& = {Too}, and P = 3 (7};) respectively. Other quantities, such as the symmetry energy, 
can be obtained from the energy density or the pressure. 

These equations are calculated for each density and proton fraction, p and y = j 
respectively, from the simultaneous solution of the field equations (obtained from the 
Euler-Lagrange ones). 



The results obtained were the following: out of 147 models only 9 of them satisfies 
all the constraints, excluding the constraint related to K x v (isospin dependence of the in- 
compressibility) [1]. These are: BSR15, BSR16, BSR17 (5% of tolerance in the PNM2 
constraint - see [1]), BSR18 (5% of tolerance in the PNM2 constraint - see [1]) [2], 
DD-F, [5] FSUGold, [6] FSUGold4, [7] FSUGZ06, [8] and TW99 [3]. Their nuclear 
matter properties are given in Table 1. It is important to mention that among these RMF 
models, 7 are nonlinear with mixing terms in the o and co fields, identified with a (*) 
and 2 are density dependent, marked with a (**). 

TABLE 1. Properties of nuclear matter at saturation density as calculated using the RMF 
models consistent with the macroscopic constraints. All entries are in MeV. 



Model 


Po 


E 


K„ 


m* 


K' 


J 


L 




BSR15* 


0.146 


-16.03 


226.82 


0.61 


512.29 


30.97 


61.79 


-252.54 


BSR16* 


0.146 


-16.05 


224.98 


0.61 


503.17 


31.24 


62.33 


-258.75 


BSR17* 


0.146 


-16.05 


221.67 


0.61 


489.45 


31.98 


67.44 


-287.31 


BSR18* 


0.146 


-16.05 


221.13 


0.61 


485.73 


32.74 


72.65 


-318.55 


DD-F** 


0.147 


-16.04 


223.32 


0.56 


758.73 


31.63 


56.00 


-285.54 


FSUGold* 


0.148 


-16.28 


229.54 


0.61 


523.93 


32.56 


60.44 


-276.07 


FSUGold4* 


0.147 


-16.40 


229.56 


0.61 


538.33 


31.40 


51.74 


-205.59 


FSUGZ06* 


0.146 


-16.05 


225.06 


0.61 


503.17 


31.18 


62.42 


-259.47 


TW99** 


0.153 


-16.25 


240.27 


0.55 


539.79 


32.77 


55.31 


-332.32 



As in the previous work [1], we see that very few models among the initial 147 
satisfy all the constraints used. We note (see Table 1) that the density-dependent models 
have values of nucleon effective mass different from the others. We also observe that 
if we include in the set of constraints the ones related to the isospin dependence of the 
incompressibility, then none of the models are consistent. Because of this, the behavior 
of these models with respect to K TW should be better investigated. 

As a final remark, we highlight the importance of the present study in the identification 
of the relativistic parametrizations which better describe the physics of both nuclear and 
pure neutron matter. 
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